Magneto-optical Kerr Effect Studies of Square Artificial Spin Ice 
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We report a magneto-optical Kerr effect (MOKE) study of the collective magnetic response of 
artificial square spin ice, a lithographically-defined array of single-domain ferromagnetic islands. We 
find that the anisotropic inter-island interactions lead to a non-monotonic angular dependence of 
the array coercive field. Comparisons with micromagnetic simulations indicate that the two per- 
pendicular sublattices exhibit distinct responses to island edge roughness, which clearly influence 
the magnetization reversal process. Furthermore, such comparisons demonstrate that disorder as- 
sociated with roughness in the island edges plays a hitherto unrecognized but essential role in the 
collective behavior of these systems. 



Arrays of lithographically fabricated single-domain 
nanoscale ferromagnets can be designed to frustrate 
inter-island magnetostatic interactions, in analogy to 
the spin-spin interactions in frustrated magnetic mate- 
rials^. A wide range of interesting behavioral can De 
observed in these artificial frustrated magnets by tun- 
ing the geometry of square^^, triangularSE', hexagonal 
(and kagomeJPplIS anc [ brickwork lattices, as well as iso- 
lated clusters^. In particular, the local moment behavior 
of these systems has shown that they are good realiza- 
tions of ice models, and a range of studies have exam- 
ined monopole excitations and drawn upon close analo- 
gies with the pyrochlore spin ice materials. 

The artificial spin ice systems, however, are subject 
to the limitations of lithography that introduce disor- 
der in the form of variation of features at the nanome- 
ter scale. While this disorder is quite different from the 
point or line defects intrinsic to an atomic lattice, it has 
the potential to influence the physics of these systems in 
interesting ways. The manifestations of disorder have 
been extensively studied in the pyrochlores and other 
frustrated magnetic materials, and off-stoichiometry and 
other structural disorder has been demonstrated to lead 
to changes in the the low temperature collective spin 
states, zero-point entropy, and other exotic phenomena 
associated with frustration HMUl Despite the demon- 
strated importance of disorder in frustrated magnetic 
materials, there has not yet been a detailed examination 
of the effects of disorder on artificial spin ice or related 
systems. 

We report magneto-optical Kerr effect (MOKE) stud- 
ies of square artificial spin ice, with in situ measure- 
ments of the global lattice magnetization. These mea- 
surements, which complement the numerous studies of 
local moments in artificial spin ice, reveal a coupling be- 
tween the collective behavior of this system and disor- 
der in the shape of the islands. We also demonstrate 
a method through which the disorder can be considered 
within micromagnetic simulations, and we provide a sim- 
ple model through which the interwoven effects of dis- 
order and magnetostatic interactions can be understood. 
The consideration of disorder opens possibilities for closer 
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FIG. 1. (color online) An SEM image of a square array at 
a — 320 nm lattice spacing with the magnetic field at an 
angle to an array axis, compared to an array of 24 ideal 
stadium-shaped islands, also at a 320 nm spacing. 



comparisons between the artificial spin ice systems and 
theoretical models, as well as for more detailed compar- 
isons with the pyrochlores, for which only the behavior 
averaged over an extensive lattice can be examined. 

Our square-ice arrays were fabricated using electron 
beam lithography, as described elsewhere 1 . The samples 
consist of 220 nm long, 80 nm wide and 25 nm thick 
permalloy (81%Ni, 19%Fe) islands with lattice constants 
ranging from 320 nm to 880 nm for the square arrays. 
The array geometry that was programmed into the elec- 
tron beam writer is shown in Figure [I] on the right. Scan- 
ning electron microscopy (SEM) of the resulting nano- 
magnet array (shown in Figure [I] on the left) revealed 
a surface roughness of ±4.3 nm on the edges of the is- 
lands, defined as the standard deviation of SEM image 
edges from an ideal island edge. As demonstrated pre- 
viously through magnetic force microscopy (MFM) stud- 
ies^, the islands are sufficiently small and elongated to 
generally behave as single-domain ferromagnets, with the 
strong shape anisotropy directing the island moments 
along their long axes. A magnetic field applied in situ 
within a MOKE magnetometer (in the longitudinal ge- 
ometrjPS) enables acquisition of full hysteresis loops dur- 
ing magnetization reversal. The 50 /am spot of an s- 
polarized HeNe laser was focused onto the arrays via an 
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optical microscope, with the sample mounted on preci- 
sion XY translation stages modified to allow sample ro- 
tation in the magnetic field. The reflected beam was 
polarization analyzed using lock-in detection to extract 
the sample-generated magnetization-dependent Kerr ro- 
tation. For an array with a lattice spacing of 320 nm 
about 40, 000 islands were simultaneously probed, while 
only ~5000 islands were probed for the largest lattice 
spacing, 880 nm. The field was incrementally ramped, 
in steps of 10 Oe near the switching fields and 100 Oe 
otherwise. The raw data were smoothed by local second- 
order polynomial regression around each point using a 
Savitzky-Golay filteM 

Figure [2] shows hysteresis loops obtained from arrays 
with 320 nm lattice spacing for different angles of the 
magnetic field with respect to one of the primary axes of 
the square lattice (as defined in Fig. [I]). The hysteresis 
loops have a strong dependence on in both shape and 
coercive field H c . For certain angles, the steepest part of 
the hysteresis curves occurs away from M = 0, such that 
M(H) near M = is only weakly sloped. Hence a defi- 
nition of H c in terms of the magnetization zero-crossing 
is problematic in this regime, and we instead define the 
coercive field as the field of maximum slope in M(H), a 
definition which enables a more consistent and unbiased 
estimate for our samples. We used a range of smooth- 
ing parameters in the Savitzky-Golay filter to estimate 
effective error bars for simulations; the error in the ex- 
perimental data was determined by the uncertainty in the 
peak fit. The extracted values of H c {0) shown in Figure [3] 
reveal a strong and non-monotonic angular anisotropy, 
with a local maximum near = 5° and a minimum at 
6 = 45°. For all lattice spacings, the qualitative features 
of the coercivity are symmetric around = 45°, as ex- 
pected due to the square lattice symmetry. We find a 
slight asymmetry in magnitude that we attribute to the 
rastering intrinsic to e-beam lithography, but the key fea- 
tures for and (90 — 0) remain qualitatively equivalent. 

We simulated the arrays' magnetic response using the 
NIST OOMMF (object oriented micromagnetics frame- 
work) code^ for an array of 24 islands in the geometry 
shown in Fig. [I] We used an OOMMF cell size of 5 x 
5x5 nm 3 , comparable to the exchange length of permal- 
1(^ 128 1 29 1 rp^ e sa t ura tj on magnetization (860 x 10 3 Am -1 ) 
and exchange constant (13 x 10 -12 Jm -1 ) are standard 
literature values^. Fig. [2] plots the resulting (M-H), 
smoothed using the same methods as for the experimen- 
tal data, and Fig. [3] displays the resulting H c (0). 

The simulation does not match the experimental data 
when using ideal stadium-shaped islands, especially the 
small minimum near = 0. This mismatch is consistent 
with previous studies showing that nominally identical 
magnetic nanostructures can have significant variations 
in their switching fields , p resumably due to shape varia- 
tions during fabrication^ 1 . Prior micromagnetic simula- 
tions have incorporated such effects by using randomly 
generated edge profiles, periodic removal of edge ele- 
ments, or edge roughness models based on experimental 
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FIG. 2. (color online) (a)- (e) Hysteresis loops comparing 
MOKE measurements (red dotted line) to simulations at 
= 0,10,15,30,45° for square arrays with 320 nm lattice 
spacing. The green dashed line and blue solid line show simu- 
lations for ideal and experimental island shapes, respectively. 
Note the significant changes with angle that are summarized 
in Figure [3] 



observation 3 - 2 -'. We incorporated island edge roughness in 
the simulation by basing the island shapes directly on 
SEM images like that of Figure [I] (the pixel size of the 
SEM images was always considerably smaller than the 
OOMMF cell size - typically around 2x2 nm 2 ). Micro- 
magnetics simulations of isolated islands show that this 
degree of shape disorder is sufficient to vary the coer- 
civity of an isolated island by ~ ±80 Oe. As shown in 
Figs. [2] and [3j the simulated magnetic response of the 



3 




10 20 30 40 

9 (Degrees) 

FIG. 3. (color online) Coercivity as a function of angle for 
a square array with 320 nm lattice spacing. Simulations were 
run on a 24-island cluster similar to that of Fig. [I] for both 
ideal and SEM-derived island shapes. Siml and Sim2 refer 
to two different simulations run for SEM images of different 
islands. 



island arrays with SEM-inferred shape disorder success- 
fully reproduces the overall shape of the experimental 
curves, including the local maximum in H c near = 5° 
and monotonic fall-off at higher angles. 

Simulations with other images or with different arrays 
of simulated islands gave qualitatively consistent results. 
We varied both the array boundaries and size of the 
arrays to test the sensitivity, and we found only a 1% 
change in coercivity when increasing the total number of 
islands in the simulation from 24 to 40 and less than 3% 
between 12 islands and 40 islands. We also observed less 
than a 1% change in coercivity at — 0° when we simu- 
lated a cluster of 24 islands chosen with different bound- 
aries from those shown in Fig. [I] The main residual dis- 
crepancy between simulations and the experimental data 
is a roughly uniform vertical shift in the simulation to 
higher overall fields. 

To elucidate the role of island-island magnetostatic in- 
teractions in the collective properties of the arrays, we 
also measured hysteresis loops for arrays with larger lat- 
tice spacings of a = 480, 880 nm, for which the inter- 
actions are much weaker than in the 320 nm array dis- 
cussed above 1 . As shown in Fig. |4^a) , the coercivity 
as a function of lattice spacing at 6 = increases with 
increasing lattice spacing for each of three different sets 
of arrays (two from one processing run and a third from 
another). Similar behavior from RefP^ (shown as array 
3) and from simulations also confirms this effect. These 
results clearly demonstrate that interactions do change 
the coercive field. As shown in Fig. [4^b), the angle- 
dependence of the coercivity at a = 880 nm shows no 
local maximum in H c near = 5°, neither in experi- 
ment nor in simulations (using the twelve-island array 
shown in the inset with SEM-derived island shapes)P^. 
Hence, we can conclude that the feature at small for 
the 320 nm lattice is associated with both island shape 
disorder and inter-island interactions and demonstrates 
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FIG. 4. (color online) (a) Variation of coercivity at 6 — as 
a function of lattice spacing for square arrays. Measurements 
on various samples (arrays 1-4) are compared with micromag- 
netic simulations (black solid squares) , comparing simulations 
of the same island outlines but with different spacing. Error 
bars for experiment are smaller than the data points, (b) 
Experiment and simulation for square arrays as a function of 
angle for the largest lattice spacing of 880 nm. The local max- 
imum in H c near — 5° that appeared for 320 nm is missing 
for the larger lattice spacing. Simulations were run on a 12- 
island array as indicated in the inset, with an appropriately 
dilated lattice spacing between the SEM islands. 



the interplay between interactions and disorder effects in 
these systems. 

We now consider how the interplay between disorder 
and the inter-island interactions lead to the increase in 
the coercive field with island separation at = and the 
local maximum in H c near = 5° for the 320 nm lattice 
constant. For small values of we can qualitatively un- 
derstand the effect of interactions among the sublattice 
of "vertical" islands, aligned with the field, by careful ex- 
amination of the micromagnetic simulations. When the 
field is swept through H C1 the magnetization of an island 
reverses suddenly when the total field (the sum of the 
external field and the field from other islands) reaches a 
critical value specific to that particular island. These mo- 
ment reversals of the vertical islands account for most of 
the change in the net magnetization of the system along 
the steep parts of the hysteresis curves. The effect of an 
island that has reversed to align with the external field 
is to enhance the external field near that island, while 
an island that has not reversed acts to reduce the mag- 
nitude of the total field acting on its neighbors. The en- 
hancement is stronger than the reduction before moment 
reversal; just before reversal, an island has complex mag- 
netization profile at the tips, whereas after reversal, its 
magnetization field is "stretched out" with strong poles. 
This asymmetry allows islands with lower intrinsic co- 
ercivity to initiate cascades of reversal^, thus decreas- 
ing the coercivity of the entire array. The cascade phe- 
nomenon, which is naturally affected by the island edge 
roughness, is clearly seen in the simulations and it ex- 
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plains at least some of the decrease in coercive field at 
small lattice spacing for all field angles. Indeed, simula- 
tions show that approximately 60% of the total decrease 
in the coercivity at = due to interactions originates 
entirely from interactions among the vertical islands. 

The sublattice of "horizontal" islands, nearly perpen- 
dicular to the applied field, plays a subtle role in altering 
the coercivity. OOMMF simulations of a 320 nm array 
of only "vertical" islands shows some flattening of the 
angle-dependence of the coercivity, relative to the coer- 
civity for the 880 nm array at small angle, but the max- 
imum is still at = 0. We therefore conclude that the 
shift of the maximum H c to a small non-zero angle, 6 max , 
is associated with the effects of the horizontal islands. A 
possible explanation is that # ma x is the angle at which 
the net effective horizontal field on the vertical islands is 
zero. According to OOMMF simulations^, at field an- 
gles of order = 5° or greater, the magnetizations of 
horizontal islands rotate in unison as the field is swept; 
near H C1 the horizontal islands are all magnetized nearly 
along their easy axes with the same orientation. This 
creates an effective horizontal field acting on the vertical 
islands that opposes and cancels the horizontal compo- 
nent of the applied field under the right conditions. 

When the external field is at zero angle, OOMMF sim- 
ulations show^ that some horizontal islands develop a 
magnetization to the left, others to the right, yet oth- 
ers go into a vortex state, a process that will be affected 
significantly by the edge profiles of the islands. In this 
situation, each vertical island experiences a different ef- 
fective horizontal field component. On average it is zero, 
but as each vertical island responds to the local field, the 
effect on H c is as though there were a small horizontal 
component to the field. Thus, a slight dip in the coer- 
civity is expected at = 0. If all the horizontal islands 
could be magnetized in the same direction while the ap- 
plied field were at 6 = 0, then the coercivity would be 
expected to drop yet more. As an experimental test, we 
temporarily applied a horizontal field to align all the hor- 



izontal islands before measurement. Subsequent MOKE 
measurements showed a decrease of H c by approximately 
10 Oe, substantially supporting the above picture. 

In summary, our data reveal the in situ collective mag- 
netization as a function of field in the 'large- array' limit 
and the influential role of island shape disorder. The re- 
sults complement previous local probe studies that have 
been very powerful in revealing short-wavelength phe- 
nomena and have demonstrated that artificial spin ice 
is a good realization of ice model physics. Given the 
high energy scales of both the island magnetic anisotropy 
and the inter-island interactions, our measurements cor- 
respond to metastable frozen states similar to those ob- 
served at low temperatures in the pyrochlore spin ice 
materials. As a result, phenomena such as the mag- 
netization steps that are observed in those materials at 
low temperatures might be expected, although a small in 
situ thermalizing fluctuation may be necessary to observe 
these effects. An interesting extension along these lines 
would be to combine MOKE techniques with dynamic 
probes such as microwave excitation of moment reorien- 
tation or local thermal excitation to above the ferromag- 
netic Curie temperature. Such a combination would add 
a new quasi-thermal aspect to the artificial frustrated 
spin ice systems and would allow a more direct compar- 
ison the pyrochlores. 
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